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Abstract: Recently we have described the discovery and complete genome sequence of a novel coronavirus 
associated with pneumonia, coronavirus HKU1 (CoV-HKU1). In this study, a detailed in silico analysis of 
the ORF lab, encoding the 7,182-amino acid replicase polyprotein in the CoV-HKU1 genome showed that 
the replicase polyprotein of CoV-HKU1 is cleaved by its papain-like proteases and 3C-like protease 
(3CL””) into 16 polypeptides homologous to the corresponding polypeptides in other coronaviruses. Sur- 
prisingly, analysis of the putative cleavage sites of the 3CL”” revealed a unique putative cleavage site. In all 
known coronaviruses, the P1 positions at the cleavage sites of the 3CL”” are occupied by glutamine. This is 
also observed in CoV-HKU1, except for one site at the junction between nsp10 (helicase) and nsp11 (mem- 
ber of exonuclease family), where the P1 position is occupied by histidine. This amino acid substitution is 
due to a single nucleotide mutation in the CoV-HKU1 genome, CAG/A to CAT. This probably represents a 
novel cleavage site because the same mutation was consistently observed in CoV-HKU1 sequences from 
multiple specimens of different patients; the P2 and P1'-P12' positions of this cleavage site are consistent 
between CoV-HKU1 and other coronaviruses; and as the helicase is one of the most conserved proteins in 
coronaviruses, cleavage between nsp10 and nsp11 should be an essential step for the generation of the 
mature functional helicase. Experiments, including purification and C-terminal amino acid sequencing of 
the CoV-HKU1 helicase and trans-cleavage assays of the CoV-HKU1 3CL” will confirm the presence of 


this novel cleavage site. 
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Coronaviruses, a genus of the family Coronaviridae, 
are large, enveloped, positive-stranded RNA viruses 
that cause respiratory, enteric, hepatic, and neurological 
diseases in humans and domestic animals. In humans, it 
has been estimated that coronaviruses [human coron- 
aviruses 229E (HCoV-229E) and OC43 (HCoV-OC43)] 
are responsible for causing about 5—30% of respiratory 
tract infections. In late 2002 and 2003, a novel clinical 
syndrome, Severe Acute Respiratory Syndrome 
(SARS), affected 30 countries on five continents with 
more than 8,000 cases and 750 deaths. A novel corona- 
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virus, the SARS coronavirus (SARS-CoV), has been 
confirmed to be the aetiological agent, and its genome 
has been completely sequenced (13, 28, 32, 36, 37, 39, 
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51-53). In early 2004, a novel human coronavirus 
associated with upper respiratory tract infections, 
human coronavirus NL63 (HCoV-NL63), was also dis- 
covered (11, 47). Phenotypically, the envelopes of all 
coronaviruses are studded with long, petal-shaped 
spikes, resulting in the appearance of a crown under the 
electron microscope. Genotypically, coronaviruses pos- 
sess the largest genomes of all RNA viruses of about 
30 kb. As a result of a unique mechanism of viral repli- 
cation, coronaviruses have a high frequency of recombi- 
nation (24, 25, 29, 31). The order of the genes in the 
genomes of all coronaviruses that encode the RNA- 
dependent RNA polymerase (pol) and the four structural 
proteins present in all coronaviruses is 5’-pol, spike (S), 
envelope (E), membrane (M), nucleocapsid (N)-3". 
Based on serological and genotypic characterization, 
coronaviruses were divided into three distinct groups, 
with HCoV-229E and HCoV-NL63 being group 1 coro- 
naviruses and HCoV-OC43 a group 2 coronavirus 
respectively (11, 26, 47). 

Recently, we have described the discovery of a group 
2 novel coronavirus associated with pneumonia, coro- 
navirus HKU1 (CoV-HKU1), and its complete genome 
sequence (48-50). In the genome of CoV-HKU1, more 
than two thirds of the genome was made up of a single 
open reading frame (ORF), ORFlab, which encodes the 
putative replicase polyprotein. In other coronaviruses, 
this polyprotein, after cleavage by proteases encoded by 
this ORF, gives rise to more than 10 proteins, of which 
many are important enzymes essential for the survival of 
the coronaviruses. In this article, we describe a detailed 
in silico analysis of this ORF in CoV-HKU1, which 
reveals a unique putative cleavage site of the 3C-like 
protease (3CL*®) in coronavirus HKU1. Other similari- 
ties and differences of this ORF in CoV-HKU1 as com- 
pared to other coronaviruses are also discussed. 


Materials and Methods 


Viral sequences. The predicted amino acid 
sequences of ORFlab in CoV-HKU1 were extracted 
from the CoV-HKU1 genome sequence (GenBank 
accession no. NC_006577) (48). The corresponding 
amino acid sequences of ORF lab in other coronaviruses 
were extracted from GenBank (HCoV-OC43, GenBank 
accession no. NP_937947; BCoV, GenBank accession 
no. NP_150073; MHV A59, GenBank accession no. 
P16342; HCoV-229E, GenBank accession no. 
NP_073549; IBV, GenBank accession no. NP_066134; 
and SARS-CoV, GenBank accession no. NP_828849). 
The amino acid sequences of the RNA-dependent RNA 
polymerases in hepatitis C virus, rabbit hemorrhagic 
disease virus and poliovirus | and their secondary struc- 


tures were retrieved from Protein DataBank (PDB ID 
1QUV, IKHV and 1RDR respectively). 

In silico analysis. Multiple alignment was _per- 
formed using ClustalX 1.83 (46). Protein family analy- 
sis was performed using PFAM (2). Secondary struc- 
ture prediction was performed using PROFsec (38). 
Three-dimensional modeling of 3CL’® of CoV-HKU1 
was performed using 3CL*® of SARS-CoV as the tem- 
plate (54). Manually corrected sequence alignment and 
homology modeling requirement were submitted via 
Deep View — spdbv 3.7 (14) to Swiss-model Protein 
Modeling web server (http://swissmodel.expasy.org/) 
(40). The 3D structure was displayed using MOLMOL 
2K.2 (22). Phylogenetic tree construction was per- 
formed using neighbour joining method with ClustalX 
1.83. 


Results and Discussion 


General Overview 

The replicase polyprotein (7,182 amino acids) of 
CoV-HKU1 is translated from ORFla (13,395 nt) and 
ORF 1b (8,154 nt). Similar to other coronaviruses, a 
slippery sequence (UUUAAAC), followed by 
sequences that form a putative pseudoknot structure, are 
present (Fig. 1) (5, 6). Translation will presumably 
occur by a -1 RNA-mediated ribosomal frameshift at 
the end of the slippery sequence. Therefore, instead of 
reading the transcript as UUUAAACGGG, it will be 
read as UUUAAACCGGG. In infectious bronchitis 
virus (IBV) and SARS-CoV, site directed mutagenesis 
experiments have confirmed the importance of the slip- 
pery sequence as well as the pseudoknot structure for 
the frameshift to occur (7, 45). 


Proteolytic Cleavage of Replicase Polyprotein 

Multiple alignments between the replicase polypro- 
tein of CoV-HKU1 and those of other coronaviruses 
reveal the orthologs of two types of proteases, papain- 
like proteases (PL®°) and 3CL’®. As in other corona- 
viruses and by multiple alignment and analysis of 
sequences around consensus cleavage sites, the repli- 
case polyprotein of CoV-HKU1 is cleaved by its PL’® 
and 3CL”° into 16 polypeptides (Table 1). The three N- 
terminal cleavage sites are putatively cleaved by the 
PL”*, whereas the others are putatively cleaved by the 
3CL™. 

The putative PL” cleavage sites in the replicase 
polyprotein of CoV-HKU1 follow the general rules 
determined by site directed mutagenesis experiments in 
other group 2 coronaviruses. In MHYV, site directed 
mutagenesis studies have shown that at the p28/p65 
cleavage site, a basic residue (arginine or lysine) in the 
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P5 position, arginine in the P2 position, and glycine in 
the P1 position are the critical amino acids for PL1”° 
recognition and processing (17). In CoV-HKU1, argi- 
nine, arginine and glycine are present in these three 
positions respectively (Table 2). In MHV, site directed 
mutagenesis studies have shown that at the p65/nsp1 
cleavage site, arginine in the P5 position, alanine in the 
P1 position, and glycine in the P1' position are the criti- 
cal amino acids for PL1”° recognition and processing 
(4). In CoV-HKU1, arginine, alanine and glycine are 
present in these three positions respectively (Table 2). 
In MHYV, site directed mutagenesis studies have shown 
that at the nsp1/p44 cleavage site, phenylalanine in the 
P6 position, glycine in the P2 position, and glycine in 
the P1 position are the critical amino acids for PL2?° 
recognition and processing (20). In CoV-HKU1, 
phenylalanine, glycine and glycine are present in these 
three positions respectively (Table 2). 

In contrast to the putative PL’® cleavage sites, analy- 
sis of the putative cleavage sites of the 3CL” in the 
replicase polyprotein of CoV-HKU1 revealed a unique 
putative cleavage site by the 3CL’® of CoV-HKU1. The 
3CL” of coronaviruses are named as such because their 
structures and substrate specificities resemble those of 
the 3C proteases in picornaviruses. In all known coro- 
naviruses, the P1 positions at the cleavage sites of the 
3CL” are occupied by glutamine, whereas the P1' posi- 
tions are occupied by aliphatic amino acids, including 
alanine, glycine, cysteine, asparagine and serine (3, 9, 
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Fig. 1. Model of ribosomal frameshift element in CoV-HKU1. It 
consists of a slippery sequence (UUUAAAC) and a pseudoknot 
structure made up by two stems and two loops. 


Table 1. Characteristics of the 16 putative polypeptides derived from replicase polyprotein in CoV-HKU1 


Fir: st amin Oa ci Position in replicase polyprotein 


Number of amino 


Polypepiides —Last amino acid?ton in spices polyprotein acid residues pharaclenshice 

p28 Met'—Gly” 222 Unknown 

p65 Tle*’-Ala®” 587 Unknown 

nspl Gly*°-Gly*** 2,029 ATR, AC, ADRP, 
PL1"", PL2"°, Y domain 

p44 Val’*°—GIn?™* 496 Hydrophobic domain 

nsp2 Ser?*-Gln’*” 303 3CL”" 

nsp3 Ser’*—Gin?™ 287 Hydrophobic domain 

nsp4 Ser?*-GiIn*” 92 Unknown 

nsp5 Ala**—Gin?”” 194 Unknown 

nsp6 Asn*?"—GIn*”? 110 Unknown 

nsp7 Ala*”-Gin*” 137 Unknown 

nsp8 Ser“*_Val*”! 14 Unknown (short peptide 

at the end of ORFla) 

nsp9 Ser“*_Gin?* 928 Pol 

nsp10 Ser***-His”® 603 Hel 

nsp11 Cys”®-Gin*” 521 ExoN 

nsp12 Ser®”-Gin*? 374 XendoU 

nsp13 Ala®*—Cys"” 299 2'-O-MT 


ATR, acidic tandem repeats; AC, acidic domain; ADRP, ADP-ribose 1"-phosphatase; PL1”°, papain-like 
protease 1; PL2"°, papain-like protease 2; 3CL””, 3C-like protease; Pol, RNA-dependent RNA polymerase; 
Hel, helicase; ExoN, 3’-to-5' exonuclease; XendoU, poly(U)-specific endoribonuclease and 2'-O-MT, S- 
adenosylmethionine-dependent 2’-O-ribose methyltransferase. 
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Table 2. Analysis of amino acid residues around putative cleavage sites of PL” and 3CL”™ in replicase polyprotein of 


CoV-HKU1 
Amino acid residues around putative cleavage sites 

Putative CoV-HKU1 MHV, BCoV, Putatively 
polypeptides Positions rin GR, Hh aks HCoV-O0C43 

pup) © 5 4 3 2 1 1 2 3° 4 5 (IL) cleaved by 
p28|p65 222/223 L Q Y I K P V L F G|V 1 
p65|nsp1 809|810 WwW F P C R K V N F — ce 
nsp1|p44 2838|2839 S L K VV LS N iL GIA 
p44|nsp2 33343335 S T S F L QS GI VK M — 
nsp2|nsp3 3637/3638 A GV K L Q S K T K R F _— 
nsp3|nsp4 39243925 E V S Q I Q S K L T D V — 
nsp4|nsp5 4013/4014 D S T V L Q A L Q S _ E F — 
nsp5|nsp6 4210/4211 A N A VM Q N N E LM P — 
nsp6|nsp7 4320/4321 S T I R L QA GV A T £E —_— 3CL”” 
nsp7|nsp9 4457/4458 S S V A V QS K DLN F _— 
nsp9|nsp10 5385/5386 K S A VM Q S V GA C V — 
nsp10|nsp11 59885989 T L P R LH C T T N L F QiC 
nsp11|nsp12 6509/6510 T F T T L QS L EN V I — 
nsp12|nsp13 6883/6884 F Y P K M Q A T N DW K — 


® 11’ amino acid residues in CoV-HKU1 are the same as those in MHV, BCoV and HCoV-OC43. 
The amino acid residues critical for cleavage by PL1”° in MHV as determined by site directed mutagenesis studies 


are highlighted. 


P4 P3 P2 Pl Pl’ P2’ p3’ 
CoV-HKUL 5985 P R H : 
MHV 5979 P R OO 
BCoV 5897 T R 
HCoV-0C43 5897 T K 
Putative cleavage site 
CoV-HKU1: 5988/5989 
MHV: 5982/5983 
BCoV: 5900/5901 
HCoV-0C43:5900;/5901 
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Fig. 2. Multiple alignment of the amino acid sequences around the putative cleavage site between nsp10 and 
nsp11 in CoV-HKU1, MHV, BCoV and HCoV-OC43. The conserved amino acids from positions P2 to P12’ in 


the four coronaviruses are highlighted. 


16, 23). In most cases, the P1' positions are occupied by 
alanine, glycine or serine. Occasionally, in group 2 
coronaviruses, the P1’ positions are occupied by cys- 
teine (30); and in rhinoviruses, by asparagines (3). In 
CoV-HKUI1, multiple alignment with the genome 
sequences of other coronaviruses revealed that all, 
except one, P| positions of putative cleavage sites of 
3CL” are occupied by glutamine. The exception lies at 
the junction between nsp10 (helicase) and nsp11 (mem- 
ber of exonuclease family), where the P1 position is 
occupied by histidine. This probably represents a novel 
cleavage site due to the following reasons. First, the 
sequence is genuine instead of a result of RT-PCR or 
sequencing errors because the same position has been 
amplified and sequenced four times in different clinical 


specimens of the index patient (48). Second, the same 
mutation was consistently found in the genome 
sequence of the CoV-HKU1 in the second patient with 
pneumonia (data not shown). Therefore, it is unlikely 
that it is due to an occasionally occurring mutation in 
just one patient. Third, in all other group 2 corona- 
viruses with genome sequences available (MHV, 
BCoV and HCoV-OC43), the P2, P1’, P2’, P3’, P4’, PS’, 
P6’, P7’, P8’, P9’, P10’, P11’ and P12’ positions of this 
cleavage site are all occupied by leucine/valine [shown 
to be important for 3CL”® cleavage (35)], cysteine, ser- 
ine/threonine, threonine, asparagines, leucine, phenyl- 
alanine, lysine, aspartic acid, cysteine, serine, 
lysine/arginine and serine respectively, which are also 
present at the corresponding positions in CoV-HKU1 
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(Fig. 2). Fourth, this amino acid substitution is probably 
due to a single nucleotide mutation in the CoV-HKU1 
genome, CAG (as in MHV) or CAA (as in BCoV and 
HCoV-OC43) (which encodes glutamine) to CAT 
(which encodes histidine). Finally, as the helicase is 
one of the most conserved proteins in coronaviruses, 
cleavage between nsp10 and nsp11 should be an essen- 
tial step for the generation of the mature functional heli- 
case. The present novel putative cleavage site, as well 
as some other atypical cleavage sites by 3CL’®, were 
not predicted by a program for prediction of 3CL”° 
cleavage sites (21). Experiments, including purification 
of the CoV-HKU1 helicase in its native form and C-ter- 
minal amino acid sequencing and trans-cleavage assays 
of the CoV-HKU1 3CL”® will confirm the presence of 
this novel cleavage site. 


Putative Polypeptides of the CoV-HKUI1 Replicase 
Polyprotein 

The replicase polyprotein of CoV-HKU1 is putatively 
cleaved by its PL’ and 3CL’® into 16 polypeptides 
(Table 1). Among these 16 putative polypeptides, p28, 
p65, nsp3, nsp4, nsp5, nsp6, nsp7 and nsp8 have no 
homologues with known functions by BLAST search. 

Similar to other coronaviruses, the N terminal of the 
putative nsp1 consists of an acidic domain (amino acids 
1-333 of nsp1). However, unlike other coronaviruses, 
there are 16 tandem copies of a 10-amino acid repeat 
near the C terminal of the acidic domain. The first 14 
tandem copies are perfect repeats of NDDENVVTGD, 
with four of the 10 amino acids being either glutamic 
acid or aspartic acid. The last two copies are imperfect 
repeats of NNDEEIVTGD and NDDQIVVTGD 
respectively. In picornaviruses, it has been found that 
there were three tandem copies of 24 imperfect amino 
acid repeats, named VPg1, VPg2 and VPg3 (10). How- 
ever, no specific functions have been found for these 
tandem repeats. Downstream to the acidic domain is 
the PL1’° (amino acids 334-549 of nsp1), with the 
characteristic catalytic dyad of cysteine and a down- 
stream histidine. Moreover, similar to the PL1™ of 
other coronaviruses, four conserved cysteine residues 
(Cys’”, Cys*’, Cys** and Cys*’) that may contain a 
Zn’** binding motif are present. Downstream to PL1”° is 
the X domain (amino acids 550—709 of nsp1) that con- 
tains putative ADP-ribose 1”-phosphatase (ADRP) 
activity [amino acids 575-685 of nsp1 belonging to 
Appr-1l-p processing enzyme family (Pfam accession 
no. PF01661)]. The putative ADRP activity was first 
described for this domain by Snijder et al. (43). In 
other microorganisms, such as Saccharomyces cerevisi- 
ae and other eukaryotes, ADRP and its functionally 
related enzyme cyclic nucleotide phosphodiesterase 


(CPDase), were important for tRNA processing (33). 
ADP-ribose 1",2"-cyclic phosphate (Appr>p) is pro- 
duced as a result of tRNA splicing. Appr>p is convert- 
ed to ADP-ribose 1"-phosphate (Appr-1"p) by CPDase. 
Appr-1"p is then further processed by ADRP. In other 
group 2 coronaviruses, both putative ADRP and 
CPDase (in NS2a) homologues are present (43). Inter- 
estingly, in CoV-HKUI1 and also SARS-CoV, only 
ADRP, but not CPDase [NS2a is not present in CoV- 
HKU 1 (48)], is present. Downstream to the X domain, 
there is a segment of 134 amino acids, also present in 
other coronaviruses, with no matches of known func- 
tions found by BLAST search. Downstream to this 
segment of unknown function is the putative PL2°° 
(amino acids 844-1140 of nsp1), with the characteristic 
catalytic dyad of cysteine and a downstream histidine. 
Similar to the PL2°° of other coronaviruses, four con- 
served cysteine residues (Cys, Cys’, Cys’ and 
Cys'*’) that may contain a Zn** binding motif are also 
present. At the C terminal of nsp1, similar to other 
coronaviruses, there is a Y domain with two hydropho- 
bic stretches and 11 conserved Cys/His residues present 
in the N terminal 889 amino acids (30). Ziebuhr et al. 
speculated that this Y domain may be responsible for 
anchoring nspl into membranes and binding cations 
(55). Downstream to nspl, p44 (496 amino acid 
residues), putatively cleaved from nsp1 by PL2?°, with a 
hydrophobic transmembrane domain, is present. 

Three-dimensional structure modeling of the putative 
3CL” (nsp2) of CoV-HKU1 using the three-dimension- 
al structure of SARS-CoV as the template revealed a 
similar structure for the monomer of the CoV-HKU1 
3CL™ as compared to the 3CL*® of SARS-CoV and 
HCoV-229E (Fig. 3). Three domains were present, 
domain I (8-101 amino acid residues), domain II 
(102-184 amino acid residues) and domain II 
(201-303 amino acid residues). Domains I and II con- 
sist of six strands of antiparallel § barrels and domain III 
is a globular structure with five @ helices. The catalytic 
dyad of the protease is located in a cleft between 
domains I (His*') and II (Cys'®). 

As in other coronaviruses, the putative Pol (nsp9) of 
CoV-HKUI1 consists of the N-terminal domain, the fin- 
gers subdomain, the palm subdomain and the thumb 
subdomain (Fig. 4). The N-terminal domain (residues 
1-371) of the Pol of CoV-HKU1 shares 52-87% amino 
acid identities with those of other coronaviruses but 
<31% amino acid identities with those of other RNA 
viruses. The function of the N-terminal domain is 
unclear. On the other hand, the fingers (motifs F and 
G), palm (motifs A—-E) and thumb subdomains of the 
Pol of CoV-HKU1 are homologous to those of other 
coronaviruses and some positive-stranded RNA virus- 
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Fig. 3. Three-dimensional structure modeling for monomer of 
CoV-HKU1 3CL”. The © helices are shown in red and f barrels 
in blue. The N and C termini are labeled N and C respectively. 
The catalytic dyad (His“ and Cys’) are indicated. 


| N terminal domain 


motif G F A B CDE 


es, such as hepatitis C virus (of family Flaviviridae) 
(1), rabbit hemorrhagic disease virus (of family Cali- 
civiridae) (34) and poliovirus 1 (of family Picornaviri- 
dae) (15) (Fig. 4). The critical amino acids, aspartic 
acids in motif A and XDD sequence in motif C, which 
have been demonstrated to be essential for the function 
of Pol in other coronaviruses and positive-stranded 
RNA viruses, are also present in Pol of CoV-HKU1. 
Similar to the helicase of SARS-CoV, the putative 
helicase (nsp10) of CoV-HKU1, which unwinds duplex 
RNA by nucleoside triphosphate hydrolysis activity, 
possesses six conserved motifs (18) (Fig. 5). Similar to 
other coronaviruses, based on the presence of conserved 
motifs, the putative helicase of CoV-HKU1 should 
belong to superfamily 1 of the three superfamilies of 
RNA helicases (12, 19). Furthermore, all nine cysteine 
residues and three histidine residues that are conserved 
in the N termini of the helicases in nidoviruses are pres- 
ent in the putative helicase of CoV-HKU1 (Fig. 5). It 
has been suggested that this cysteine/histidine region 
constitutes a Zn’*-binding domain, which controls the 
activities of the catalytic domain (42). Similar to the 
helicases of other nidoviruses and most other helicases 
of superfamily 1, the helicase of CoV-HKUI1 probably 
unwinds the RNA duplex in a 5‘ to 3’ fashion (41, 44). 
The putative nspl11 possesses a putative 3’-to-5’ 
exonuclease (ExoN) domain of the DEDD superfamily 
(56). The putative nsp12 possesses a putative poly(U)- 
specific endoribonuclease (XendoU) domain (27). The 
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Fig. 4. Schematic representation of the predicted secondary structures of the RNA-dependent RNA polymerases in 
CoV-HKU1, HCoV-OC43, HCoV-229E and IBV and the secondary structures of the homologous parts of the RNA- 
dependent RNA polymerases in hepatitis C virus, rabbit hemorrhagic disease virus and poliovirus 1 (large solid bar, o- 
helix; small empty bar, B-sheet). In the multiple alignment of conserved motifs, amino acid residues identical in the 
four coronaviruses are marked by stars and conserved amino acid residues in all seven viruses are shaded. 
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Fig. 5. Multiple alignment of amino acid sequences of the helicases of CoV-HKUI, MHV, BCoV, HCoV-OC43, 
SARS-CoV, HCoV-229E and IBV. The nine cysteine residues and three histidine residues that are conserved in the N 
termini of the helicases in nidoviruses are highlighted. The six conserved motifs are also depicted. 
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putative nsp13 possesses a putative S-adenosylmethio- 
nine-dependent ribose 2'-O-methyltransferase (2°-O-MT) 
domain of the RrmJ family (8). These three putative 
enzymes, as well as ADRP and CPDase, are enzymes in 
RNA processing pathways. ExoN, XendoU and 2'-O- 
MT are enzymes in a small nucleolar RNA processing 
and utilization pathway, in contrast to the pre-tRNA 
splicing pathway that ADRP and CPDase belong to. At 
the moment, no experiments have been performed to 
confirm the activities of these putative RNA processing 
enzymes. Further studies are required to elucidate the 
exact roles of these putative viral enzymes in the corre- 
sponding viruses. 
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